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ABSTRACT 

We present evidence that the WASP-14 exoplanetary system has misaligned orbital and stellar- 
rotational axes, with an angle A =— 33.1°± 7.4° between their sky projections. The evidence is based 
on spectroscopic observations of the Rossiter-McLaughlin effect as well as new photometric observa- 
tions. WASP-14 is now the third system known to have a significant spin-orbit misalignment, and all 
three systems have "super- Jupiter" planets (Mp > 3 Mj up ) and eccentric orbits. This finding sug- 
gests that the migration and subsequent orbital evolution of massive, eccentric exoplanets is somehow 
different from that of less massive close-in Jupiters, the majority of which have well-aligned orbits. 
Subject headings: stars: individual (WASP-14) — planetary systems: individual (WASP-14b) — 
techniques: spectroscopic — techniques: photometric 



1. INTRODUCTION 

Close-in giant planets are thought to have formed at 
distances of several AU and then migr ated inward to 
their current locations (|Lin et al.l [l996) . The mecha- 
nism responsible for the inward migration of exoplan- 
ets is still debated. Some clues about the migration 
process may come from constraints on the stellar obliq- 
uity: the angle between the stellar spin axis and the or- 
bital axis. The sky projection of this angle, A, can be 
measured by observing and interpreting the anomalous 
Doppler shift during the tran sit of a planet, known as the 
Rossiter-McLaughlin effe c t (iMcLaughlml 119241: iRossiter 
19241: iQueloz et all l200Cj _ Winn et alj 120051: lOhta et al 
20051: Gaudi & Winn 2007). Some of the proposed mi- 



gration pathways would produce large misalignments (at 
least occasionally) while others would preserve the pre- 
sumably close alignment that characterizes the initial 
condition of planet formation. 

For example, theories that invoke migration of the 
planet through interactions with the gaseous protoplan- 
etary disk predict small spin-orbit angles, and that ini- 
tial spin-orbi t misalignments and eccentricities should be 
damped out (iLin et ail 1 19961: IMoor head fc Adamsll2008t 
iLubow fc Qgilviel [2001). On the other hand, impul- 
sive processes such as c lose encounters betw een planets 
()Chatteriee et all 120081: iFord fc R,asioll2008ft or dynami- 
cal relaxation (jJuric fc Tremaing|2008[ ) should drive sys- 
tems out of alignment. The Kozai mechanism also pro- 
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duces large orbital tilts (Wu & Murray 2003, Fabrycky & 
Tremaine 2007). Ultimately the hope is that the predic- 
tions of migration theories ca n be compared with an e n- 
semble of measurements of A (|Fabrvckv fc Wirmll2009f ). 

In this paper we add the transiting exoplanet WASP- 
14 b to the growing collection of systems for which the 
projected spin-orbit angle has been measured. WASP- 
14 is a relatively bright (V — 9.75) F5V star which was 
discovered by the Wide- Angle Search for Planets (Super- 
WASP) to unde rgo periodic trans its by a Jovian planet 
every 2.2 days (| Joshi et alj 120091 hereafter J09). The 
planet is among the most massive of the known tran- 
siting exoplanets, with Mp — 7.3 Mj up , and it has a 
measurably eccentric orbit (e = 0.090 ± 0.003), which 
is unusual among the hot Jupiters. J09 also reported 
a measurement of the spin-orbit angle, A = — 14_^ de- 
grees, which is consistent with zero, but also allows for 
the possibility of a significant misalignment. In the fol- 
lowing section we describe our spectroscopic and photo- 
metric observations of WASP-14, made in an attempt to 
refine the measurement of A. In §[3] we present evidence 
for a large spin-orbit misalignment based on our radial- 
velocity measurements obtained during transit. We sum- 
marize the results of our joint analysis of our photometric 
and spectroscopic monitoring in § |4j and present tenta- 
tive evidence of an emerging trend between spin-orbit 
misalignment, and the physical and orbital characteris- 
tics of close-in exoplanets. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Radial Velocity Measurements 

We observed the transit predicted by J09 to occur on 
2009 J une 17 using the H igh-Dispersion Spectrometer 
(HDS, iNoguchi et all 120021 ) on the Subaru 8.2 m Tele- 
scope atop Mauna Kea in Hawaii. We obtained spectra 
of WASP-14 through an iodine cell using the I2b spec- 
trometer setting and a 0'.'8 slit, providing a resolution 
of approximately 60, 000. We started our observing se- 
quence just after evening twilight, about 20 min before 
the predicted time of ingress. We continued our observa- 
tions until 2.5 hr after egress when the star set below 20° 
elevation. For most of our observations we used exposure 
times of 5 min, yielding a signal-to-noise ratio (SNR) of 
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Fig. 1. — Relative radial velocity measurements of WASP- 14 
as a function of orbital phase, expressed in days since midtransit. 
The symbols are as follows: Subaru (circles), Keck (squares), Joshi 
et al. 2009 (triangles). The lower panel shows the residuals after 
subtracting the best-fitting model including both the Keplerian 
radial velocity and the Rossiter-McLaughlin effect. 

100-120 pixel -1 at 5500 A, the central wavelength of the 
range with plentiful iodine absorption lines. At high air- 
mass we increased our exposure times to 10 min. 

We also obtained 8 radial velocity measurements of the 
G2V star HD 127334 on the same night, 2009 June 17. 
HD 127334 is a long-term target of the California Planet 
Search. Keck/HIRES radial velocity measurements over 
the past 3 years show that the star is stable with an rms 
scatter of 2.5 m s" 1 . With HDS we used exposure times 
ranging from 60 to 120 seconds, resulting in SNR ranging 
from 110-120 pixel" 1 at 5500 A. 

We also obtained out-of-transit (OOT) radial velocities 
of WASP-14 using the High-Resolution (HIRES) spec- 
trometer on the Keck I telescope starting in July 2008. 
We set up the HIRES spectrometer in the same man- 
ner that has been used consiste ntly for the California 
Planet Search (|Howard et al.ll2009f ). Specifically, we em- 
ployed the red cross-disperser and used the I2 absorption 
cell to calibr ate the instrumental re sponse and the wave- 
length scale (jMarcv fc Butlerlll992f ). The slit width was 
set by the 0'.'86 B5 decker, and the typical exposure times 
ranged from 3-10 min, giving a resolution of about 60,000 
and a SNR of 140-250 pixel" 1 at 5500 A. 

For the spectra obtained at both telescopes, we per- 
formed the Dopp ler analysis with the algorithm of 
iButler et ail (|1996h . as updated over the years. A clear, 
Pyrex cell containing iodine gas is placed in front of the 
spectrometer entrance slit. The dense forest of molecu- 
lar lines imprinted on each stellar spectrum provides a 
measure of the wavelength scale at the time of the obser- 
yation, as well as the sh ape of the instrumental response 
(|Marcv fc Butlerlll992]) . The Doppler shifts were mea- 
sured with respect to a "template" spectrum based on 
a higher-resolution Keck/HIRES observation from which 
the spectrometer instrumental response was removed, as 
far as possible, through deconvolution. We estimated the 
measurement error in the Doppler shift derived from a 
given spectrum based on the weighted standard devia- 
tion of the mean among the solutions for individual 2 A 
spectral segments. The typical measurement error was 



Fig. 2. — Top Panel: Relative photometry of WASP-14 during 
the transit of 2009 May 12. Bottom Panel: Residuals from the 
best-fitting transit light curve model. 

1.0-1.7 m s" 1 for the Keck data and 6 m s" 1 for the Sub- 
aru data. The RV data are given in Table [T] and plotted 
in Figure [T] 

2.2. Photometric Measurements 

We observed the photometric transit of 2009 May 12 
with the University of Hawaii 2.2 m (UH2.2m) telescope 
on Mauna Kea. We used the Orthogonal Parallel Trans- 
fer Imaging Camera (OPTIC), which is equipped with 
two Lincoln Lab s CCID128 ortho gonal transfer array 
(OTA) detectors (jTonrv et al.lll997t ). Each OTA detector 
has 2048x4096 pixels and a scale of 0'.'135 pixel" 1 . OTA 
devices can shift accumulated charge in two dimensions 
during an exposure. We took advantage of this charge- 
shifting capability to create large square-shaped point 
spread functions (PSFs) that permit longer exposures be- 
fore reaching saturation dHowell et "aT1l2003l ; iTonrv et al.1 
[20051: 1 Johnson et al.ll2009fl . 

We observed the transit of WASP-14 continuously for 
5.5 hr spanning the transit. We observed through a cus- 
tom bandpass filter centered at 850 nm with a 40 nm 
full-width at half- maximum. We shifted the accumulated 
charge every 50 ms to trace out a square-shaped region 
25 pixels on a side. Exposure times were 50 s, and were 
separated by a gap of 29 s to allow for readout and re- 
freshing of the detectors. Bias subtraction and flat-field 
calibrations were applied u s ing cu stom IDL procedures 
described by I Johnson et al.1 (|2009f ). 

The only suitable comparison star that fell within the 
OPTIC field of view is a V = 12.1 star ~ 6 arcmin- 
utes to the Northeast. The fluxes from the target and 
the single comparison star were measured by summing 
the counts within a square aperture of 64 pixels on a 
side. Most of the light, including the scattered-light halo, 
was encompassed by the aperture. We estimated the 
background from the outlier-rejected mean of the counts 
from four rectangula r regions flanking each of the stars 
(| Johnson et al.|[2009f ). As a first order correction for vari- 
ations in sky transparency, we divided the the flux of 
WASP-14 by the flux of the comparison star. The tran- 
sit light curve is shown in Figure [H and the photometric 
measurements and times of observations (HJD) are listed 



3 



in Table H 



3. DATA ANALYSIS 
3.1. Updated Ephemeris 



The first step in our analysis was to refine the esti- 
mate of the orbital period using the midtransit time de- 
rived from our OPTIC light curve. We fitted a transit 
m odel to the light curve based on the analytic formulas 
of iMandel fc A gol (2002) for a quadratic limb-darkening 
law. The adjustable parameters were the midtransit time 
T t , the scaled stellar radius i?*/a (where a is the semima- 
jor axis), the planet-star radius ratio Rp/ R+, the orbital 
inclination i, the limb darkening coefficients u\ and u-i , 
and two parameters k and mo describing the correction 
for differential airmass extinction. The airmass correc- 
tion was given by 



mobs + mo + kz 



(1) 



where m bs is the observed instrumental magnitude, z is 
the airmass, and m cor is the cor rected magnitude that is 
compared to the transit model (jWinn et al.l l2009). 

We used the rms of the OOT measurements as an 
estimate for the individual measurement uncertainties. 
We did not find evidence for significant tim e-correlated 
noise using the time-averaging method of iPont et ahl 
(200(| . We fitted the light curve model and esti- 
mated our parameter uncertainties using a Markov Chain 
Monte Carlo algorithm (MCMC; iTegmark etldl 12004 
lFordll2005t [Gregorvll200"5h . The results for the midtran- 
sit time, and the refined orbital period (using the new 
midtransit time and the midtransit time given by J09) 
are 

T c = 2454889.8921 ±0.00025 (2) 
P = 2.2437704 ±0.0000028 days. (3) 

The other derived lightcurve parameters were consistent 
with those reported by J09. 

3.2. Evidence for Spin-Orbit Misalignment: A Simple 

Analysis 

Figure [2] shows our Subaru/HDS and Keck/HIRES RV 
measurements made near transit, after subtracting the 
best-fitting Keplerian orbital model. The RVs measured 
just after ingress are redshifted with respect to the Ke- 
plerian orbital velocity. We interpret this "anomalous" 
redshift as being due to the blockage by the planet of 
the blueshifted limb of the rotating stellar surface. We 
therefore conclude that the planet's orbit is prograde. 

In addition, the anomalous redshift persists until about 
1 hr after midtransit. This is evidence for a misalignment 
between the orbital axis and stellar rotation axis. Were 
A = 0, the midpoint of the transit chord would be on 
the projection of the stellar rotation axis, and therefore 
the anomalous Doppler shift would vanish at midtransit, 
in contradiction of the data. Thus we can conclude that 
the orbit of WASP-14b is inclined with respect to the 
projected stellar spin axis. In the next section we make 
a quantitative assessment of A. 

7 We allowed both coefficients to be free parameters, subject to 
the conditions < «i + U2 < 1 and «i > 0. 



3.3. Global Analysis of Radial Velocities and 
Photometry 

We simultaneously fitted a parametric model to our 
OPTIC light curve and to the the 4 sets of RV 
data: Subaru/HDS and Keck/HIRES (this work), and 
OHP/SOPHIE and NOT/FIES (J09). To make our anal- 
ysis of the RM effect largely independent of J09, we did 
not include the RVs gathered by J09 during transits. The 
photometric aspects of the model were given in § 13.11 The 
RV model was the sum of the radial component of the 
Keplerian orbital velocity, and the anomalous velocity 
due to the Rossiter-McLaughlin effect. To compute the 
latter, we used the "RM calibration" procedure of Winn 
et al. (2005): we simulated spectra exhibiting the RM 
effect at various orbital phases 8 , and then measured the 
anomalous radial velocity AVr of the simulated spectra 
using the same algorithm used on the actual data. We 
found the results to be consistent with the formula 



1.124- 0.395 



3.5 km s 



> (4) 



where 5f is the instantaneous fractional loss of light dur- 
ing the transit and v p is the radial velocity of the occulted 
portion of the stellar disk. 

The 18 model parameters can be divided into 3 groups. 
First are the parameters of the spectroscopic orbit: the 
period P, the midtransit time T t , the radial- velocity 
semiamplitude K, the eccentricity e, the argument of 
pericenter cu, and velocity offsets for each of the 4 differ- 
ent groups of RV data. Second are the photometric pa- 
rameters: the planet-to-star radius ratio Rp/R*, the or- 
bital inclination i, the scaled stellar radius i?*/a (where a 
is the semimajor axis), the 2 limb-darkening coefficients, 
and the out-of-transit flux and differential extinction co- 
efficient. Third are the parameters relevant to the RM 
effect: the projected stellar rotation rate vsini* and the 
angle A between the sky projections of the orbital axis 
and the stellar rotation axis [for illustrations of the ge- 
ometry, see Ohta et al. (2005), Gaudi & Winn (2007), or 
Fabrycky & Winn (2009)]. 

The fitting statistic was 




P - 2.243770 d 



0.0000028 d 



where /j(obs) are the relative flux data from the OPTIC 
light curve and afj is the out-of-transit rms. Likewise 
Vj (obs) and a v j are the RV measurements and uncertain- 
ties. For a v j we used the quadrature sum of the measure- 
ment error and a "jitter" term of 4.4 m s -1 , which was 
taken from the empirical calibration of Wright (2004). 

8 For the template spectrum, which should be similar to that 
of WASP-14 but with slower rotation, we used a Keck/HIRES 
spectrum of HD3681 {T eff = 6220 K, [Fe/H]= +0.08; vsinu = 
2.8 ± 0.5 km s" 1 IValenti fc Fischerll207j5h . 
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Fig. 3. — Relative radial velocity measurements made during transits of WASP-14. The symbols are as follows: Subaru (circles), 
Keck (squares), Joshi et al. 2009 (triangles). Top Panel: The Keplerian radial velocity has been subtracted, to isolate the Rossiter- 
McLaughlin effect. The predicted times of ingress, midtransit and egress are indicated by vertical dotted lines. Middle Panel: The residuals 
after subtracting the best-fitting model including both the Keplerian radial velocity and the RM effect. Bottom Panel: Subaru/HDS 
measurements of the standard star HD 127334 made on the same night as the WASP-14 transit. 



The final term enforces the constraint on the orbital pe- 
riod based on the new ephemcris described in the previ- 
ous section. 

As before, we solved for the model parameters and 
uncertainties using a Markov Chain Monte Carlo algo- 
rithm. We used a chain length of 5 x 10 6 steps and 
adjusted the perturbation size to yield an acceptance 
rate of ~40%. The posterior probability distributions 
for each parameter were approximately Gaussian, so we 
adopt the median as the "best-fit" value and the stan- 
dard deviation as the l-er error. For the joint model fit 
the minimum \ 2 1S 291.4 with 295 degrees of freedom, 
giving xt — 0.99. The contributions to the minimum \ 2 
from the flux data and the RV data were 246.1 and 45.3, 
respectively. The relatively low value of the RV contri- 
bution compared to the number of RV data points (64) 
suggests that 4.4 m s _1 is an overestimate of the jitter 
for this star, and that consequently our parameter errors 
may be slightly overestimated, but to be conservative we 
give the results assuming a jitter of 4.4 m s . 

For the main parameter of interest, the projected spin- 
orbit angle, our analysis gives A = 1 ± 7.4° (Figure [3J . 
Thus, the WASP-14 planetary system is prograde and 
misaligned, as anticipated in the qualitative discussion of 
§ 13.21 Our measurement of A agrees with the value mea- 
sured by J09 (— 14^4 deg), but with improved precision 
that allows us to exclude A = with high confidence. 



We find the projected stellar rotational velocity to be 
vsini* = 2.80 ± 0.57 km s _1 . This value is somewhat 
lower than, but consistent with, the values determined 
by J09 from line broadening wsin = 3.0 ± 1.5 km s _1 , 
from their RM analysis vsini* = 4.7 ± 1.5 km s , and 
from our SME analysis fsinz* = 3.5 ± 0.5 km s" 1 . This 
agreement among the rotation rates provides a consis- 
tency check on our analysis. The best-fitting parameters 
and their uncertainties are listed in Table [3] 

We found no evidence for another planet or star in the 
WASP-14 system. To derive quantitative constraints on 
the properties of any distant planets, we added a single 
new parameter 7 to our model, representing a constant 
radial acceleration. A third body with mass M3 <C M*, 
orbital distance 83 > a and inclination would produce 
a typical radial accleration 



and our result is 7 = 2.0 ±1.4 cm s 1 d , or 1.01 ± 
0.72 Afj up (5 AU)- 2 . 

4. SUMMARY AND DISCUSSION 

We present new photometric and spectroscopic mea- 
surements of the WASP-14 transiting exoplanetary sys- 
tem. By combining a new transit light curve, several 
Keck/HIRES RV measurements made outside of tran- 
sit, and most importantly, Subaru/HDS RVs spanning a 
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Fig. 4. — The spin-orbit configuration of the WASP-14 planetary 
system. The star has a unit radius and the relative size of the planet 
and impact parameter are taken from the best-fitting transit model. 
The sky-projected angle between the stellar spin axis (diagonal 
dashed line) and the planet's orbit normal (vertical dashed line) is 
denoted by A, which in this diagram is measured counter-clockwise 
from the orbit normal. Our best-fitting A is negative. The 68.3% 
confidence interval for A is traced on either side of the stellar spin 
axis and denoted by a\. 



transit, we have measured and interpreted the RM effect. 
By modeling the RM anomaly we find that the projected 
stellar spin axis and the planetary orbit normal are mis- 
aligned, with A = -33.1°± 7.4°. 

Of the 13 transiting systems with measured spin- 
orbit angles, only 3 have clear indications of spin-orbit 
misalignments. The other two cases besides WASP-14 
are XO-3 ( Hebrard et al. 2008, Winn et al. 2009) and 
HD 80606 (|Gillonl 120091: IPont et al.ll200l ; Winn et al. 
2009c in prep). It is striking that all 3 tilted systems 
involve planets several times more massive than Jupiter 
that are on eccentric orbits, and that none of the sys- 
tems with eccentricities consistent with circular or with 
masses smaller than 1 Mj up show evidence for misalign- 
ments (Figure [S]). 

In addition to the three known super- Jupiters with 
inclined orbits, there are also two eccentric, mas- 
sive exoplanets with small proje c ted spin-orbit a ngles : 
HP 17156 b (iFischer et all 120071 IB arbieri et all [20071 : 
ICochran et alj|2008t iNarita et alJ [ 2009h and HD147506 
(Bak os et all 120071: IWinn et all l2007t lLoeillet et all 
|2008| ). However, neither of these cases presents as strong 
an exception to the pattern as it may seem. The mea- 
surement of A in both cases was hampered by the poor 
constraint on the transit impact parameter, which causes 
a strong degeneracy between A and usini (Gaudi & 
Winn 2007). It should also be kept in mind that the 
measured quantity A is only the sky-projected spin-orbit 



angle, and that the true angle of one or both of those 
systems may have a stellar rotation axis that is inclined 
by a larger angle along our line of sight. 

It was already known that the orbits of massive planets 
are systematically diffe rent from the orbits o f less massive 
planets. For example, Wr ight et al. (2009), building on 
a previous finding by iMarcv et al.1 (|2005l ). showed that 
planets with minimum masses Mp sin i > 1 Mj up typ- 
ically have lower orbital eccentricities than those with 
minimum masses smaller than 1 Mj up . While sub- 
Jovian-mass planets have eccentricities that peak near 
e = with a sharp decline toward e = 0.4, those with 
Mp sin i > 1 Mj up have eccentricities that are uniformly 
distributed between e = and 0.55. 

The tendency for misaligned orbits to be found among 
massive planets on eccentric orbits does not yet have a 
clear interpretation. It may seem natural for inclined or- 
bits and eccentric orbits to go together, since both incli- 
nations and eccentricities can be excited by few-body dy- 
namical interactions, whether through the Koz ai effect, 
(|Fabrvckv fc Tremaind 120071: IWu et all I2007D planet- 
planet scattering ( Juric &; Tremaind 12008 . see, e.g.,), or 
scenarios combining both of these phenomena (Nagasawa 
et al. 2008). However, the mass dependence of these and 
other mechanisms for altering planetary orbits needs to 
be clarified before any comparisons can be made to the 
data. 

The misalignment of the WASP-14 planetary system, 
along with the previously discovered misaligned systems, 
have offered a tantalizing hint of an emerging trend 
among the orbital and physical properties of close-in, 
transiting exoplanets. However, trends seen in small data 
sets can often be misleading. To bring this picture into 
better focus, a more sophistic ated analysis of the extant 
data, following the example of iFabrvckv fc Winnl (2009) 
may be better. And as with all astrophysical trends, 
observations of a larger sample of objects will provide 
a much clearer picture than any statistical analysis of 
a smaller sample. Thus, additional RM observations of 
transiting systems are warranted, with particular atten- 
tion paid to trends with orbital eccentricity and planet 
mass. 
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TABLE 1 

Radial Velocity Measurements of WASP- 14 



Heliocentric Julian Date (HJD) RV Uncertainty Tclcscope a 

(m s _1 ) (m s _1 ) 



2454667 8(1421 


-139.4 


1.0 


K 


2454672.81824 


-1008.4 


1.3 


K 


2454673 83349 


955.3 


1.3 


K 


2454999 76227 


151.6 


6.3 


s 


2454999 76665 


138.8 


5.2 


s 


2454999 77D91 


133.6 


5.9 


s 


2454999 77517 


112.0 


5.4 


s 


2454999 77943 


110.2 


5.5 


s 


2454999 7929D 


97.9 


5.8 


s 


2454999 7971 6 


89.8 


5.4 


s 


2454999.80142 


67.8 


7.1 


s 


2454999 806011 


65.8 


6.0 


s 


2454999 81 026 


54.3 


5.4 


s 


2454999 81 452 


49.0 


5.2 


s 


2454999 81878 


34.6 


5.3 


s 


2454999 82901 


4.9 


4.7 


s 


2454999 83327 


-6.9 


4.5 


s 


2454999 84471 


-38.3 


5.4 


s 


2454999 84898 


-46.0 


4.8 


s 


2454999 85771 


-70.3 


4.9 


s 


2454999 861 97 


-76.8 


5.1 


s 


2454999 86623 


-92.5 


4.7 


s 


2454999 87049 


-111.3 


4.5 


s 


2454999 87477 


-114.5 


5.7 


s 


2454999 87904 


-124.9 


6.0 


s 


2454999 88330 


-131.0 


5.8 


s 


2454999 88757 


-148.0 


5.5 


s 


2454999.89183 


-161.0 


6.0 


s 


2454999 89610 


-154.8 


5.6 


s 


2454999.90037 


-175.6 


5.8 


s 


2454999.90464 


-187.6 


5.4 


s 


2454999.90890 


-190.4 


6.7 


s 


2454999.91317 


-215.8 


7.1 


s 


2454999.92129 


-228.0 


8.0 


s 


2454999.96315 


-330.4 


8.8 


s 


2454999.96743 


-331.1 


7.8 


s 


2454999.97170 


-350.1 


9.2 


s 


2454999.98117 


-372.5 


5.6 


s 


2454999.98891 


-378.0 


5.4 


s 


2454999.99665 


-398.5 


5.9 


s 


2455000.00825 


-416.6 


5.7 


s 


2455000.01598 


-429.6 


6.4 


s 


2455014.86287 


965.8 


1.4 


K 


2455015.91393 


-812.7 


1.7 


K 



a K: HIRES, Keck I 10m telescope, Mauna Kea, Hawaii. S: HDS, Subaru 
8m telescope, Mauna Kea, Hawaii. 



TABLE 2 
Relative Photometry for WASP-14 



Heliocentric Julian Date (HJD) Relative Flux 



2454963.85021 


1.00064 


2454963.85113 


1.00127 


2454963.85204 


1.00024 


2454963.85296 


1.00086 


2454963.85387 


1.00115 


2454963.85478 


0.99986 


2454963.85569 


1.00183 


2454963.85660 


1.00005 



Note. — The full version of this table is avail- 
able in the online edition, or by request to the au- 
thors. 



TABLE 3 
System Parameters of WASP-14 



Parameter Value 



Orbital Parameters 

Orbital period, P [days] 2.2437704 ± 0.0000028 

Mid-transit time, Tt [HJD] 2454963.93676 ± 0.00025 

Velocity scmiamplitude, [m s _1 ] 989.9 ±2.1 

Argument of pcriccntcr, ui [degrees] 253.10 ±0.80 

Orbital eccentricity, e 0.0903 ± 0.0027 

Velocity offset, 7 FIES [m s" 1 ] -4989.5 ± 3.4 

Velocity offset, 7sophie [m s" 1 ] -4990.1 ± 3.0 

Velocity offset, 7hires [ m s_1 ] 107.1 ± 2.1 

Velocity offset, 7hds [m s _1 ] 7.7 ± 2.5 



Spin-orbit Parameters 

Projected spin-orbit angle A [degrees] 

Projected stellar rotation rate v sinu [km s -1 ] 



-33.1° ± 7.4° 
2.80 ± 0.57 



